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A discontinuous generalization of the standard map, which arises naturally 
as the dynamics of a periodically kicked particle in a one dimensional infi- 
^ ' nite square well potential, is examined. Existence of competing length scales, 

namely the width of the well and the wavelength of the external field, in- 
troduce novel dynamical behaviour. Deterministic chaos induced diffusion is 
Q | observed for weak field strengths as the length scales do not match. This is 

related to an abrupt breakdown of rotationally invariant curves and in partic- 
■ ular KAM tori. An approximate stability theory is derived wherein the usual 

standard map is a point of "bifurcation". 
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O 1 I. INTRODUCTION 

IT) 

o 

The construction and study of area preserving mappings has led to a deeper understand- 
ing of apparently complex dynamics, especially of Hamiltonian chaos. The maps range from 
abstract models |1J such as the cat maps and the baker map to more generic ones. One 
of the most well studied of such generic mappings is the standard map [0, which has also 
been investigated extensively in its quantum version The classical map on the cylinder 
displays a range of dynamical behaviours, and when completely chaotic, diffusive random 
walks in momentum take place. Recent experiments using trapped ultra-cold sodium atoms 
in pulsed laser fields have probed this model, as a kicked rotator, and verified the central 
phenomenon of quantum localization of momentum diffusion ||. 

In this Letter we study the classical dynamics of a simple generalization of the standard 
map that naturally arises from the dynamics of a particle trapped in a one dimensional in- 
finitely deep well. The virtue of this generalization is that on the introduction of a competing 
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length scale novel dynamical behaviour manifests, leading to large scale diffusive processes 
even for small external field strengths. This places the standard map within a family of 
generally discontinuous area preserving mappings; the points where there is continuity cor- 
responding to the usual standard map. 

We study the novel dynamical changes by simple methods exploiting separation of time 
scales in non-autonomous equations. The failure of the Poincare-Birkhoff theorem leads to 
the formation of entirely stable islands and cantori or entirely unstable and chaotic orbits. 
It is also to be pointed out that this can be one of the methods of controlling or enhancing 
chaos. 

The classical motion of a particle in an infinite square well potential in the presence of a 
uniform monochromatic external field has been studied previously || . In this case diffusion 
was observed over a limited range of energies and there is no issue of competing length 
scales. Also, no analytic form of the stroboscopic mapping is derivable. We have considered 
below a space variation of the external field that is also pulsed in time, the pulsing leading 
to analytically derivable mappings. More recently a special case of this system was studied 
H and quantized, where (in the chaotic regime) an important effect was found, namely 
derealization of eigenstates as opposed to the well known exponentially localized states 
of the kicked rotator. Our independent and parallel studies are a generalization and in 
this Letter we present the classical aspects as we feel that these deserve a more complete 
understanding. 

Recent advances in semiconductor physics, makes it possible to construct wells on atomic 
scales. Also, developments in investigating the quantum nature of electrons in a finite quan- 
tum well have been achieved, using tunnel-current spectroscopy, when the corresponding 
classical system shows chaotic behaviour (see references therein). We hope that the 
study of novel dynamical behaviours in simple models, as the generalization of the standard 
map, will reflect on some aspects of such experimentally realizable systems. 



II. KICKED PARTICLE IN A WELL 

We consider a particle of mass m, trapped in an one dimensional infinite square well 
potential V sq (x) of width 2a. There is an external field V(x) that is periodically pulsed with 
period T. The Hamiltonian is 

H = + V sq (x) + V(x)±J (j - |) (1) 

where x,p are the position and momentum of the particle respectively. We consider below 
V(x) = e cos(27rs/A). Here e and A are the amplitude and wavelength of the external field 
respectively. The effect of the pulsed field in H is the same as that of an infinite number of 
travelling waves with identical amplitudes and frequencies which are multiples of 2tt/T (the 
pulse frequency). 

The stroboscopic map, relating the dynamical variables immediately after successive 
kicks, can be derived in a standard manner and results in the following dimensionless form: 

X n+1 = (-l) M "(X„ + P n ) + (-l) M " +1 Sgn(P n )M n 

P n+ i = (-l) M "P n + (A/27r)sin(2vr J RX n+1 ). 
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where M n 



Sgn(P n )(X n + P n ) + ± 



, is the number of bounces of the particle at the walls 
during the interval between the nth kick and the (n + l)th kick. Here Sgn(..) and [..] stand 
for sign and integer part of the argument respectively. The dimensionless quantities are 
defined as 

X.* P. = 2£, K — ~~"~~T~; (3) 
2a 2am am\ A 

while the dynamical variables x n , p n are the position and momentum of the particle just 
after the nth kick. Although Sgn(P n ) is discontinuous and undefined for P n = 0, one can 
see easily from the map that either of the values (±1) can be taken for Sgn(P n ) as this does 
not alter the dynamics. We refer to the mapping (Q) as the well map. Note that K and R 
are the only two effective parameters of the well map. Here K is related to the strength of 
the external field; R is the ratio of the width of the well to the wavelength of the external 
field. We do not impose any constraints on R, for example those while considering standing 
waves in a cavity, but allow it all values. Evidently, \X n \ < 1/2 or \x n \ < a as the particle 
motion is confined between two rigid walls. 

We now relate the well map to a generalized standard map. The motion of a free particle 
inside the potential well V sq (x) is related to the motion of a free rotator, the difference being 
one of boundary conditions. Both free motions are same unless the particle hits the walls. 
The correspondence between them is made explicit in the following way. Let us consider the 
dynamics of a free rotator in discrete unit time steps, which is the map relating successive 
angles 6 n and angular momenta J n . 

n +i =9n + J n (mod 1) ; J n+1 = J n . (4) 

Here the motion is confined to a cylinder [—1/2, 1/2) x (—00, 00). Similarly the discretized 
free motion in a well of unit width is 

X n+1 = {-l) Mn {X n + P n ) + (-l) M " +1 Sgn(P n )M n ; P n+1 = (-l) M "P n (5) 

Let us denote S n = (9 n , J n ), and W n = (X n , P n ) . If we have the same initial condition for 
both the maps ([|) and (|5]) i.e., Wo = So, we have the following relation after the nth iterate: 

W n = (-l) M S n (6) 

where M = X)?=o -^i> is the total number of bounces of the particle after (n — 1) iterates. In 
the rest of the paper we consider the mapping 



Jn+i = J n + {K/2ir)sm{2irR0 r , 

On+l = 9n + Jn+1 (mod 1). 



(7) 



It is easy to see that we have the same relation @ between the well map @ and the map 
defined in ([!]) (more precisely the time reversal of this map). 

We refer to the map defined in ([?]) as the Generalized Standard Map (GSM). From the 
above description it is clear that quantitative dynamical features are same for the well map 
and the GSM. In attempting to analyze the dynamics of the well map, it is sufficient to 
understand the dynamics of the GSM. Also the map (0) is easier to handle than the well 
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map, as it is only a slightly generalized version of the well studied standard map (R = 1) [||. 
When it is stated that the standard map is a one parameter system the implicit assumption 
is that there is only one length scale. Here there are naturally two length scales whose 
ratio is the cause of many interesting effects. Both the maps (fj) and (0) have reflection 
symmetry about their origin. Below we discuss various dynamical behaviour of the GSM in 
more detail. Throughout this paper we consider R to assume positive real values only. 



III. DYNAMICAL FEATURES 

For the standard map (GSM with R = 1), when K -C 1, the dynamics is quasi-regular 
and the phase space is filled with a large number of KAM tori which are rotationally invariant 
circles || extending across the phase plane. These closed loops C encircle the cylinder and are 
invariant, meaning that TC = C, where T is an area preserving continuous transformation 
on the cylinder (for instance Eqn.(|7]) with R — 1). KAM tori are the principal barriers for 
the unstable orbits to diffuse in the phase space. As K increases there is a smooth transition 
from regular to chaotic behaviour with a reduction in the number of KAM tori. Dynamical 
changes of KAM tori with increase of K are discussed in [[|. At K ~ 1 all the KAM tori 
disappear from the phase space [HI; the dynamics is chaotic and diffusive for K 3> 1. Note 



that the GSM is periodic in J and 9 with unit period. For integer values of R, it is continuous 
and takes the form of standard map with K —>■ KR. In this case dynamical transitions are 
identical with that of the standard map. As the effective perturbation here is KR, higher 
integer values of R also make the dynamics more chaotic. 

When R assumes non-integer values the situation is entirely different from the earlier 
case, particularly when K is small. Shown in Fig. [I] are typical phase space structures of the 
GSM for R ^ 1 and < R < 2 at small K. It is seen that no KAM tori exist even when R is 
slightly away from unity. In general, no KAM tori appear in the phase space when R departs 
from integer values. This is of course due to the discontinuity in the map when R assumes 
non-integer values. In other words, there are no closed loops C encircling the cylinder such 
that T'C = C, where T is the transformation defined by GSM with non-integer values of R. 
As can be seen easily the discontinuity in T arises at 9 = —1/2 (or 1/2). 

This discontinuous dynamics implies the abrupt breakup of KAM tori and a failure of 
the Poincare- Birkhoff scenario. The modifications that may occur are little understood, and 
Chirikov in [11 points out that the piecewise linear sawtooth map (p n +i = p n +K x n , x n+ i = 



%n + Pn+i mod 1) when A' is a non-integer such that —4 < K < is one such system where 
extremely complicated locally stable motions occur and "it is not at all clear what could be 
a meaningful description, if any, of this apparently trivial model". The particle in a well 
naturally gives rise to a non-linear generalization of the sawtooth map and below we show 
how a simple method may give us significant local stability information that explains the 
dramatic differences observed in Fig. [1|. 

In addition Fig. [I] distinguishes the dynamics between the cases R ^ 1 and R ^ 1. 
The former case is more chaotic and the phase space is filled with unstable orbits, while 
many regular and chaotic regions are seen in the latter. Thus the entire standard map may 
be regarded as being poised at a point of bifurcation when regarded as a function of the 
parameter R. However, it is seen that when K is large, there is no such apparent difference 
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in dynamics between the above two cases since the phase space becomes chaotic and diffusive 
in both the cases. 



A. Stability of the tori 

In the absence of the general KAM theorem for irrational tori and the Poincare-Birkhoff 
theorem for rational tori, we expect to apply some kind of local analysis to study deviations 
from integer values of R. Indeed we will see that the integer values of R could be thought of 
as bifurcation points for entire orbits. We noted earlier that even for small values of K no 
KAM tori exist in the phase space of the GSM when R departs from unity. To explain the 
stability changes we look at the behaviour of a single irrational KAM torus, which exists at 
R = 1, as R varies. A typical observation (Fig. |2|) shows that KAM torus becomes unstable 
and chaotic for R & 1 while it is replaced by a chain of stable islands for R ^ 1. 

We begin by performing a linear stability analysis for a KAM tori with respect to R. We 
assume below that K is small enough for such an orbit to exist. Let us denote a GSM (with 
R = k, positive integer) KAM orbit as {8 n , J n }. This satisfies the following mapping: 

Jn+l = Jn + (K/2n) sin(27rfc0 n ) 

(8) 

9 n+1 = n + J n +i (mod 1). 

Also consider an orbit {8 n , J n } of the GSM with R = k + fi and <C 1, with the same 
initial conditions as the above KAM orbit. Let us introduce the differences A6 n = 9 n — 6 n ; 
AJ n = J n — J n . Then it is possible to write an exact mapping equation for A9, AJ as 

AJ n+1 = AJ n + (K/2rr) {sin (2tt(A: + //)(0 n + A0 n )) - sin(27r£;0 n )} 

[i>) 

A9 n+ i = A6 n + AJ n+ i. 

Expanding in terms of A6 n and retaining first order terms leads to a time dependent force 
and the non-autonomous linear set of equations: 

AJ n+1 = AJ n + K(k + /i) cos(27r(A; + jj)6 n )A6 n + A n 
A6 n+1 = A9 n + AJ n+ i 



with 



A n = — fsin(2vr(A; + n)0 n ) - sin(27r£;0 n )} . 



2ti 

The behaviour of such linear non-autonomous equations can be quite complex (compare 
for instance the Mathieu differential equation which also arises in linear stability analysis). 
For our analysis on the stability of the tori we wish be intuitive and derive rough but useful 
estimates. Although we never expand R about an integer, here we assume that this excur- 
sion is small, so that we can treat the non-autonomous stability equations perturbatively. 
Appealing to the method of averaging we simply replace the time dependent force by its 
time average. The limitation of this is pointed out further ahead. 
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Also since the motion on the KAM torus is ergodic in 9, we replace the time average by 
an equivalent space average with uniform measure. Thus 



coa(2irRB n ) » g(R) = cos(2ttRx) dx (11) 

7-1/2 



This procedure leads us to: 

AJ n+1 = AJ n + (K/ir)sin(7r(k + n))A9 n + A n 

(12) 

A9 n+1 = A9 n + AJ n+1 . 

Eqn. fll^D is a simple linear map and its stability can be seen from the corresponding Jacobian 
matrix which gives the stability condition | 2 + (K/tt) sin(7r(A; + //)) | < 2. This leads to the 
following window of stability: 

-% < (-l) k sm(nfi) < 0. (13) 
K 

The above stability window implies that if k is odd KAM tori which exist at R = k are 
stable for /i ^ and unstable for ji ^ 0; converse is the case if k is even. Although \i is 
small in the above analysis, our numerical observations satisfy the stability condition ( |T3"D 
even for large values of \x (see Fig. |2|, which is with k — 1). This explains large phase space 
regions in Fig. |I|. 

Before dwelling more in this region of phase space we divert attention to those where 
this procedure has failed. The chaotic regions around ( J = 0, 9 = 0) in Fig. [I] for instance 
presents such a case. In the interior of the original 0/1 resonance (at R = 1) even a reversal 
of the stability criteria is observed. The hyperbolic fixed point at R = 1 persists for R ^ 1 
and is the genesis of the chaos. The smooth stable and unstable manifolds of the hyperbolic 
fixed point that exist for R = 1 and small K cannot exist for R ^ 1 as they would then 
represent rotationally invariant curves. Therefore the presence of the hyperbolic points 
necessarily implies chaos. 

Within the framework of the stability analysis done earlier what has failed around the 
hyperbolic fixed points is the assumption of uniform measure. If we think of the simple 
pendulum we are moving from energetic rotational motions (KAM) to slow oscillations of 
large amplitude and clearly the time is spent preferentially around the turning points. Also 
the replacement with time averages will fail as the time scales involving fast 9 motions and 
slow A J motions become comparable. The latter effect becomes important as we move into 
the 0/1 resonance region. We note that the question of the existence of two different time 
scales is initial value dependent and is present in the standard map as well, i.e., it does not 
result of the presence of the additional parameter R. 

The modified measure can be approximated well by considering the simple pendulum. 
We now take k = 1 for simplicity and consider only excursions of R from unity. The term 
cos(27ri?^ n ) in Eq. (|T0| ) is replaced with the ergodic average: 

cos{27rRd n ) « g{R,9 c ) = / 1 ; / / (14) 

Jdc J E — cos(27rx) Je c J E — cos(27rx) 
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Here — 1 < E < 1 is the scaled energy and 9 C is the turning point given by cos(27r# c ) = E. 
This is shown in Fig. ^ and the point where g crosses zero from above ought to be a point 
where stability is recovered. 

This then explains the uniform chaos around the hyperbolic fixed point for R ^ 1. This 
also explains the stability around the wall {9 — 1/2, J — 0) in all cases. However the 
interior of the resonance is not accessible to the simple theory. Remarkably Fig. |3] predicts 
a recovery of stability at 8 C ~ 0.15 which is seen when R = 0.95 while stability is recovered 
at a much higher value of 6 C for R = 1.05 and is not seen in Fig. |3|. This fact we attribute 
to the failure of the averaging procedure in analyzing non-autonomous equations and have 
no other recourse than Eq. (|H]) itself. 

As we have noted, hyperbolic fixed points generate chaos when R ^ 1. The lack of 
chaos in the "meat" of the phase space for R ^ 1 where there were KAM tori are therefore 
regions where there is a total absence of hyperbolic points. The Poincare-Birkhoff theorem 
concerning the breakup of rational tori into an equal number of hyperbolic and elliptic 
periodic points is clearly violated. In fact the modified scenario seems to be the creation of 
only elliptic fixed points for R ^ 1 and only hyperbolic points for R ^ 1. To illustrate this 
we have shown in Fig. |] the fate of points on the lines J = 1/3 and J = (\^E — l)/2. While 
in the former case three elliptic islands emerge prominently, in the latter one can see chains 
of elliptic islands with the clear hierarchy of the number of islands 5, 8, 13, 21, 34 . . . deriving 
from the Fibonacci sequences generating closer approximations to the golden mean. 

These island chains alternate as the ratios approach the golden mean. As a limiting 
case we may think of a golden mean chain, but we have of course not rigorously proved its 
existence. There appear to be two types of orbits in the stable regions: those that form 
an island chain and those that meander in the interstitial spaces between island chains on 
presumably a fractal set. These would then be examples of strange non-chaotic sets in 
Hamiltonian mechanics. More investigations into these simple systems seems warranted. 

We have thus acquired a rather detailed understanding, using simple methods, of the 
dramatic stability changes that accompanies the introduction of an additional length scale 
or discontinuity. We may also speculate that this may prove a good way of controlling or 
enhancing chaos. While our study has been for Hamiltonian systems dissipative systems 
may also display such a behaviour. 

B. Lyapunov Exponents 

The Lyapunov exponents can be calculated using the Jacobian matrix of the GSM, which 
at nth iterate we denote as M n . The stability changes induced by R ought to reflect on the 
exponent somewhat like a phase transition effect on an order parameter. Shown in Fig. ^ is 
a typical calculation of the positive Lyapunov exponent A + for K = O.ln. The orbit chosen 
is a KAM torus for R = 1 and we, from our earlier discussion, expect it to be unstable for 
R < 1 and stable for R > 1 but not too large. It is seen that A + ~ for R ^ 1 in the 
range shown, which corresponds to the regular motion due to the stable islands formed from 
the orbit; for R < 1, A + > which indicates the chaotic behaviour of the orbit. Also the 
orbit has a maximum value of the Lyapunov exponent at R ~ 0.5, and is hence maximally 
unstable at this point. We now explore this analytically. 

The Jacobian has trace 
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Tr(M n )| = \2 + KRcos(2nRe n )\ 



(15) 



and hence |Tr(M n )| > 2 for R < 1/2 as \9 n \ < 1/2. We remind the readers that this range 
of 9 n is the entire configuration space and 9 n = ±1/2 represent the wall boundaries. Since 
the Jacobian is only a function of the angle at a given phase space point (and not also the 
action), we conclude that the Jacobian has real eigenvalues throughout the phase space. In 
other words, the system is hyperbolic for R < 1/2. This implies that there are contracting 
and expanding real directions or alternatively stable and unstable manifolds throughout 
phase space. It is worth remarking that very few dynamical models are known exactly to be 
completely chaotic or hyperbolic. These include the sawtooth maps ]T2|, the baker maps, 
flows on surfaces of constant negative curvature. The standard map even for large values 
of the parameter K is not proven to be completely hyperbolic. Thus in this context, we 
place the fact that GSM has a parameter range for which it is completely hyperbolic for all 
positive values of the parameter K. 

The similarity of the GSM, when R < 1/2, with the piecewise linear sawtooth maps 
is made clear by the following linear approximation. We approximate the force f(6 n ) = 
sm(2TrR8 n ) in the GSM by a monotonic linear function 2 sin(7ri?)6 l n . This leads to |Tr(M„) | « 
1 2 + {K/tx) sin(7T.R)|. The Lyapunov exponents in this approximation are given by 

A± = ln(r ± Vr 2 - 1) (16) 

where r = 1 + (K/27r) sin(7ri?). It is clear from Fig. |5| that the approximated exponent fits 
fairly well. Similar agreement was seen for a wide range of initial conditions and different 
values of K as well. This shows the validity of the above linear approximation for f(8 n ) when 
R < 1/2, which enables us to understand the gross behaviour of the Lyapunov exponents 
in this regime. We note that the usual linear approximation f(9 n ) ~ 2nR6 n near the origin 
was found not to be as good as the above approximation. 

Although this study has been essentially that of another area-preserving map, we feel 
that these novel effects are worth emphasizing, especially since the standard map has been 
realized in the laboratory in the form of kicked cold sodium atoms, and the transitions to 
chaos in the quantum regime have been studied as noted above. Implementation of the 
competing length scales in this experiment could lead to novel observable transport effects. 
Theoretical studies of the quantum dynamics of this system are currently underway. 
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Figure Captions 

Fig 1 : Phase space of the GSM (with R ^ 1) at K = O.Itc. No KAM tori are seen in any 
of the cases. 

Fig 2 : Shown is an orbit with initial condition So = (0,0.34), at K = O.ln, which corre- 
sponds to a KAM torus when R — 1. For R ^ 1, the orbit becomes unstable and chaotic, 
while for R 1 it becomes chain of stable islands. 

Fig 3 : The function g(R, 8 C ) which when positive determines the instability of some GSM 
orbits. 

Fig 4 : Orbits with initial conditions on the line J = 1/3 (top) and J = (V5 — 1)/2 (bottom) 
for K = 0.1 and R = 1.05. 

Fig 5 : Positive Lyapunov exponent (dots) of the orbit shown in Fig. |2|. The exponent is 
positive for R < 1 with maximum at R ~ 0.5 and zero for 1 < R < 2. Solid curve is the 
approximation (|16|) for R < 1/2. 
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